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Abstract. The particle-scale dynamics of granular materials have commonly 
been characterized by the self-diffusion coefficient D. However, this measure 
discards the collective and topological information known to be an important 
characteristic of particle trajectories in dense systems. Direct measurement of 
the entanglement of particle space-time trajectories can be obtained via the 
topological braid entropy S'braid? which has previously been used to quantify 
mixing efficiency in fluid systems. Here, we investigate the utility of S'braid 
in characterizing the dynamics of a dense, driven granular material at packing 
densities near the static jamming point ^j. From particle trajectories measured 
within a two-dimensional granular material, we typically observe that S'braid is 
well-deflned and extensive. However, for systems where (j) > 0.79, we flnd that 
S'braid (like D) is not well-defined, signifying that these systems are not ergodic 
on the experimental timescale. Both iSbraid cind D decrease with either increasing 
packing density or confining pressure, independent of the applied boundary 
condition. The related braiding factor provides a means to identify multi-particle 
phenomena such as collective rearrangements. We discuss possible uses for this 
measure in characterizing granular systems. 
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1. Introduction 

Recently, there have been extensive efforts to draw parahels between dense granular 
systems, foams, emulsions, and glassy molecular systems [l][2]. For such systems 
approaching the glass/jamming transition, the dynamics slow down due to an 
increasing intrication of the available phase space, with rearrangements taking the 
form of cage jumps ^jjs]. In this regime, particles are trapped by their neighbors at 
short timescales and resume diffusive behavior at long timescales due to eventual loss of 
correlation in the succession of jumps. Such jumps involve cooperatively-rearranging 
regions which induce dynamical heterogeneities [9, 10 , and the particles pass each 



other only rarely. The growth in size of these rearranging regions is believed to be 
associated with the rigidity and loss of ergodicity exhibited by fragile glasses. 

Unlike glassy molecular systems, granular materials permit the tracking of 
individual particles. The average dynamics have traditionally been characterized by 
the mean squared displacement, cr^(r). When the experimental timescale is long 
enough compared to the viscous timescale, permitting a complete decorrelation in the 
particle motion, it is possible to calculate the self-diffusion constant D. This coefficient 
is a 5m^/e-particle measure in which all information regarding relative motions of 
particles is discarded. In contrast, recent experiments and simulations [3|[8| pT||T8] 
have shown that the dynamics near the glass transition become heterogeneous in 
both space and time. As such, it is necessary to take relative motions into account 
when formulating the correct average description of the state. To this end, multi- 
particle scalars have been developed, including the four-point susceptibility X4 |12f[14] . 
However, to our knowledge, no scalars have been introduced which characterize mixing 
effects, which may be relevant to quantifying the cooperative motion of particles. 

The central observation of this paper is that the trajectories of a two-dimensional 
assembly of grains take the form of a braid when plotted in a space-time diagram, as 
shown in Figure [l] Each strand of the braid represents the space-time trajectory 
of a single particle in the system, with the length of the braid corresponding to 
time and the entanglement of the strands arising from dynamics. The topological 
braid entropy, S'braid^ is a measure of the degree of this entanglement. This approach 
has been successfully used to analyze the mixing efficiency in two-dimensional fluid 
flows |T£-24 . Specifically, the degree of lagrangian chaos created by the motion of 



rods in a fluid is quantitatively related to the topology of the corresponding braid. 
We anticipate that this notion might also be applicable to characterizing caging and 
dynamical heterogeneity in granular systems, by providing a quantitative measure of 
how close a system is to departing from ergodicity. 

Here, we experimentally investigate the behavior of the average topological braid 
entropy extracted from the long-time dynamics of particles in a dense, driven granular 
material, examining the dependence on boundary conditions (constant pressure, CP, 
and constant volume, CV) and the inter-particle friction coefficient ji. We analyze 
the dynamics with two methods: a single-particle method (the self-diffusion constant, 
D) and a topological multi-particle method (the braid entropy, 6'braid)- We find that 
both D and iSbraid are sensitive to the packing density (j) and pressure P of the system 
in the same qualitative way. However, at high (j) and P, both D and 5'braid are 
inaccessible due to insufficiently long experimental timescales or ill-definedness. In 
contrast, the braiding factor can be readily computed and offers an instantaneous 
view of the magnitude and intermittency of the rearrangements in the system. 
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2. Experiment 

We conduct experiments on a two-dimensional granular assembly supported upon a 
nearly-frictionless horizontal air table and agitated by bumpers around the boundary. 
The granular material is bi-disperse, consisting of large and small particles with 
diameter of rL = 83.6 mm and rs = 55.8 mm and mass of tul = 8.1 g and ms = 3.5 g, 
respectively. The concentration of small/large grains is fixed at Ns/Nl = 2, so as to 
occupy similar areas. The apparatus is shown in Figure [l] 

The boundary condition and the inter-particle friction ja can be changed for each 
experimental run. The boundary condition is established by the configuration of the 
piston (mpiston = 95 ms). By attaching weights to the piston via a low-friction 
pulley, the system is confined under constant pressure (CP) with a range of pressures 
5.7 < P < 80mN/m. The constant volume (CV) boundary condition is established 
by fixing the piston to the table, where a range of packing densities = 0.72 to 0.81 
is explored by removing units of particles (two small and one large), maintaining the 
relative concentration. The inter-particle friction is selected by changing the material 
around the particles' outer edge, with /ii =0.1 for PTFE (Teflon) wrapping, /i2 = 0.5 
for bare polystyrene, and /is = 0.85 for rubber. 

Around the boundary, an array of bumpers agitates the system on three sides. 
Bumpers are triggered pairwise via a computer at a high frequency (/ = 10 Hz) 
to keep particles in motion despite energy dissipated in collisions. A triggered 
pair consists of a bumper and its corresponding bumper on the opposite wall. We 
have previously measured 25 that the effect of this driving system is to provide a 
thermal-like bath which maintains the granular material at constant average kinetic 




Figure 1. (Color online) (a) Schematic and (b) photograph of apparatus. The 
piston provides constant pressure via pulleys and weights, or constant volume by 
fixing its position to the surface of the table, (c) Entanglement of space-time 
trajectories for = 20 particles. 
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Figure 2. Representative diffusion relations cr'^(r) for constant pressure (CP) 
experiments with /i2 and (j) = 0.74 to 0.80. The dashed hue represents a slope of 
unity. The diffusion coefficient, D, is calculated between the two + symbols on 
each curve. 



energy. More information on the driving and kinetics of the system can be found in 
references 25 - 27 . 

Images are taken Tcamera = 2 or 5 seconds apart depending on the (/)-dependent 
dynamical timescales [26]. We record the positions of particles and map trajectories in 
time using unique particle identifiers 26 . The unique identifiers allow a relatively large 
^camera, SO wc cau Confidently track the trajectories of individual particles and thereby 
be certain that topological dynamics do not erroneously result from the exchange of 
particle identities due to tracking mistakes. The packing density <p is calculated as 
(/) = {7rr^/Vi)i with the particle radius and Vi its local Voronoi cell area, calculated 
with Voro++ [28]. The average is taken over all particles in the central 20% area of 



the table to reduce ordering effects induced by the boundaries 25 29 



3. Results 



3.1. Self- diffusion 

The diffusion behavior is quantified using the mean-squared displacement 

a\T) = {Ax,{tj + + Ayi{tj + r)^),,, (1) 

where the average (•) is performed over all particles i and all times tj. As can be seen 
in Figure [2] the slope of cr^(r) is not constant for all timescales r. It is useful to define 
a local curve cr^(r) oc r", where the value of a characterizes the average individual 
dynamics associated with that r. For short timescales (r < td)^ the trajectories 
remain sub-diffusive, with a < 1. This arises because the particles are locally 
confined, or caged, by neighboring particles [3j[4j. At longer timescales (r > r^)). 
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each particle has experienced multiple cage-breaking events and thus resumes diffusive 
behavior. For timescales at which a = 1, we can measure the two-dimensional self- 
diffusion coefficient D = |^cr^ ('?")• Finally, if the diffusion process is fast enough, a 
typical particle covers a distance of many particle diameters within the duration of the 
experiment. In such cases, the exponent a can again fall below unity at long times, 
due to the finite size of the experiment. This effect is visible in Figure [2] for the lowest 
(j) at which we perform experiments. 

In general, it is possible to calculate D using cr^(r) for experiments which are 
of long enough duration to reach the a = 1 regime. For this system, we find that 
this regime is reached even at values of (j) near static random loose packing [26]. A 
convenient means to measure D is to specify a tolerance on a which defines a diffusive 
regime, marked by + symbols in Figure [2j We calculate D by fitting to the largest 
continuous region where the smoothed derivative of cr^(r) has a slope \a — 1| < 0.1. 
The diffusive timescale td is the timescale at which the particle motion first achieves 
diffusive behavior. With increasing packing density ^, td becomes longer, as seen in 
Figure [2] For systems with large (j) and P, the range of r over which (r) was used 
to find D (marked with crosses in Figure |2| becomes less than a decade. 

Given the constraint that |a — 1| < 0.1 and the finite duration of our experiment 
~ 10^ s, it is not always possible to calculate a diffusion constant for systems with 
large P or 0. One such example (0 = 0.80, P = 80mN/m) is shown in Figure [2] A 
more conservative tolerance specification on a would further decrease the experiments 
for which we could calculate a meaningful value of D. In many experiments where we 
would wish to quantify the self-diffusion of particles, D is in fact poorly-defined. 



3.2. Braid entropy 



The braid entropy is a measure of the degree of entanglement of trajectories 2Qp2p3 



To define it, first consider an imaginary rubber band enclosing the strings at the 
bottom of the braid shown in Figure As we slide the rubber band upwards, its 
length L{t) will tend to grow, as it is caught on the strings and is not allowed to 
traverse them. If the trajectories are chaotic, then L{t) will grow exponentially; the 
growth rate is then related to the Lyapunov exponent of a dynamical system. The 
braid entropy iSbraid is defined as the asymptotic growth rate of 

5'braid = lim ^logL(t), (2) 

t^oo at 

maximized over the choice of rubber band. This rubber band viewpoint is well- 
illustrated by taffy pullers and dough-kneading devices, where material is made to 



stretch exponentially by repeated folding 30 . The length L{t) is called the braiding 
factor. 

When trajectories arise from a continuous dynamical system, then the braid 
entropy is a lower bound on the topological entropy of the system. The topological 
entropy is a measure of the loss of information about the identity of trajectories — it 
is an entropy in the sense of information theory. As more trajectories are included. 



the braid entropy converges to the topological entropy 31 . However, in a granular 
medium such as the one considered here, there is no underlying continuous dynamical 
system: the braid entropy is not necessarily related to some intrinsic topological 
entropy. 

We follow the method described in 23 to calculate a good approximation to the 



braid entropy. This is a symbolic approach based on earlier work ^ 33 where a 
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Figure 3. The braiding factor logL(t) for the same experimental runs shown in 
Figure [2] with (f) = 0.74 to 0.79. The inset shows log L(t) divided by its value at 
the maximum time, T. 

braid is first converted into its algebraic group representation for a particular choice 
of projection axis, and then used to measure the growth of hypothetical loops similar 
to the rubber band described above. We calculate the braiding factor L{t) using the 
trajectories of = 20 particles. Because particles can enter and leave the region 
captured by the camera, this choice of number provides a compromise between having 
too many particles with missing points in the trajectory, and too few particles with 
the maximum trajectory duration. We find that this choice of N is sufficient for L{t) 
to capture the dynamics of the whole system, as will be discussed in more detail below. 

The algebraic group representation is constructed from a list of the crossing times 
and identities of the particles involved in each crossing. Linear interpolation is used to 
determine the instant of crossing along the projection line, and thus obtain the braid 
generators. The experimental data is sufficiently well-resolved for this method 23 to 
accurately capture crossings. 

The braid entropy 5'braid (see Equation [2| is measured from the growth rate of 
the logarithm of the braiding factor. To obtain a single value for a particular run, we 
calculate a linear fit to logL(t). This entropy will be well-defined if logL(t) meets 
three conditions: grows linearly in time, is an extensive quantity, and is independent 
of the choice of projection axis. We can directly check each of these conditions. 

In Figure [sj we show log L{t) for systems identical to Figure |2j omitting the 
{(j) = 0.80, P = 80 mN/m) system as too few persistent crossings occur. We find 
log L{t) grows, on average, linearly with time for (j) < 0.79. In the inset, log L{t) is 
divided by its value at the maximum time, to highlights how the degree of linearity 
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Figure 4. The S'braid per particle, shown for 122 and CP systems, with the axis 
of projection ^ = 0. The arrow shows increasing cf) of the system. 




Figure 5. A plot of logL(t) for two systems with /i2 where the axis of projection 
is rotated between = [0,7r/2], where (a) is (f) = 0.77, P = 12 mN/m and (b) is 
(j) = 0.78, P = 18 mN/m. 



present at different values of (j). 

To test the extensivity of 6'braid, we repeat the calculation described above for 
different numbers of particle trajectories, and test whether S/N is a constant. Figure[4] 
shows the results of this test for a series of CP experiments at different 0. For A/" > 15, 
we observe that S/N is quite flat, independent of boundary condition, /i, ^, and P. 
This observation is consistent with iSbraid being an extensive quantity. For very small 
ensembles of trajectories, < 15 (less than 10% of the system size), we find that the 
braid is too sparse to be representative of the dynamics of the system. In addition, it 
is difficult to test for extensivity in systems with low (j). In this regime, particles have 
a larger diffusion constant and are therefore more likely to enter and leave the region 
being monitored by the camera. Therefore, it is not possible to find N ^ 25 particles 
which all satisfy the constraint on the minimum duration of continuous trajectories 
required to make a calculation of 5'braid- 

We compute L{t) for 9 projections along axes oriented from = to 7r/2, and 
compare the results in Figure [5] for two characteristic systems, both confined at CP and 
with /i2 but at different (j) and P. For the example at low ^, we observe that log L{t) is 
quite linear, even on short timescales, and that the choice of projection axis has little 
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effect on the slope of log L{t) or 5'braid- However, in the higher-^ example, log L{t) 
grows more quickly for the {0 = 7r/2)-projection than for the {0 = 0)-projection, where 
the difference in entropy is ~ 2%. Increasing N decreases the difference in entropy 
between projections. In the limit of long experiment duration, such differences in 
calculated entropies would be expected to vanish. 

In examining the runs at large (j) and large P, we note the occurrence of 
two important features in log L{t): fiat regions and steep increases. In Figure [s}) 
(representative of such systems), we observe that log L{t) may remain constant for an 
extended period of time, as in this case for the interval 0.7 < t < 1.2 x 10"^ s. Such a flat 
L{t) curve signifies that particles are not exchanging neighbors. The second feature 
to note in Figure|5]3 is the sharp increases in log L{t) beginning around t ^ 1.5 x 10^ s, 
for all choices of projection axis. Such steps correspond to a large number of crossings, 
no matter the chosen orientation, and therefore signify a collective rearrangement of 
particles (cage-breaking event). Another such event can be seen in Fig. [s] just after 
t = 1 X 10^ sec. Due to these events, the slope of log L{t) can be much larger on 
short timescales. As such, the assumption that log L{t) grows linearly in time is not 
valid for short timescales. However, we find Shraid is well-defined for systems where 
(0 < 0.79, P < 50) as log L{t) is linear on the experimental timescale, extensive, and 
independent of axis of projection. 



3. 3. Comparison 

Based on prior work [4l|26,34 , we expect D to decrease as (f) increases. Indeed, in 
Figure [6^b, D decreases with increasing (j) and P, independent of boundary condition 
(CP and CV). Similarly, we find that 5'braid decreases with increasing (j) and P, as 
shown in Figure [HJ^d, respectively. Error bars represent the average standard error of 
'S'braid from the set of 9 projections on axes oriented from 9 — to 7r/2, where the 
magnitude of the standard error of 6'braid is 10% and does not depend on and 
P. We observe little dependence on the inter-particle friction coefficient /j. on D and 
-Sbraid^ though qualitatively D decreases with increasing ja. 

As discussed earlier, an experimental measurement of D using requires a 
constraint on the value of a to be used to find which was selected to be |q^ — 1| < 0.1. 
For systems with high P, this constraint was not met and therefore D is poorly defined, 
as shown in Figure |6^b. Similarly, 5'braid is also poorly-defined for high P systems 
as log L(t) no longer grows linearly within the experimental timescale. For systems 
with (j) > 0.79 and P > 50 mN/m, neither D nor ^braid are useful for measuring the 
dynamics even with long experimental timescales of 5 x lO^s (~ 14 hr). 



4. Discussion & Conclusion 



We have investigated the behavior of the global braid entropy 5'braid associated with 
the long-time dynamics of a dense granular assembly, and find it to be a well- 
defined quantity. Notably, 5'braid was observed to be extensive, making it a promising 
candidate to play a role in the elaboration of a proper equation of state for dense 
granular systems. The breakdown in the measurement of ^braid coincides with the 
loss of diffusive dynamics, which is to be expected due to the association of the onset 
of caging behavior with a loss of ergodicity. The measurement of both ^braid and 
the single-particle diffusion coefficient D provide a quantitative measure of the degree 
to which slower individual dynamics also results in poorer mixing. In spite of the 
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Figure 6. Measured values of D as a function of cf) and P in (a) and (b), and 
likewise jShraid as a function of cf) and P in (c) and (d), respectfully. Error bars 
denote the magnitude of the standard error. The black x on the axis signifies 
systems for which D or S'braid is poorly defined for that particular cf) or P. (e) 
A logarithmic plot of S'^raid D for each system, where the arrow denotes 
increasing cf). 



fact that D is di single-particle measurement, we find that it is a good predictor of 
the degree of mixing. One benefit of characterizing the average dynamics with ^braid 
instead of D is that no arbitrary fitting parameters are necessary: the method is fully 
specified by the trajectories themselves. 

In addition to providing a measure of the average behavior, the braiding factor 
L{t) provides a simple way to identify temporal events governing the dynamics (the 
crossings). In particular, L{t) displays increasing intermittency as the packing density 
or pressure is increased, with successions of plateaus and jumps indicative of collective 
neighborhood changes, or structural rearrangements. In fact, this departure from 
linearity can provide a measure of the chaoticity/ergodicity of the dynamics. The 
duration required to reach an acceptable degree of linearity, independent of projection 
axis, is a measurement of a characteristic dynamical timescale of the system. From 
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this point of view, the experimental timescales over which we have investigated our 
dense granular systems were comparable or larger than this equilibration time, with the 
exception of the highest P and <p systems we tested. For our experimental timescale, a 
loss of ergodicity was observed for systems with close proximity to the static jamming 
transition. 

This investigation represents a first step towards a more systematic and controlled 
characterization of the braid entropy in two-dimensional granular systems. The 
braiding factor L{t) may be a means of obtaining system trajectories which could be 
analyzed with the recently developed thermodynamics of histories 35 . Quantitative 
benchmarking of 5'braid in thermal and glassy systems, together with a complete 
comparison with the usual dynamical susceptibilities are promising directions for 
future investigations. 
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